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Abstract 

There has been a recent interest in considering Quantum Field Theories in which Lorentz 
Invariance is broken in the UV sector. However attention has been mostly limited to dis- 
persive theories. In this work we provide the generalized settings for studying dissipation. 
Unitarity is preserved by coupling the original fields to additional (heavy) fields which 
induce the dissipation. Starting with Lagrangians breaking LI in the UV, we learn that 
dissipative effects unavoidably develop in the effective theory. We then covariantize these 
Lagrangians in order to address the trans-Planckian question of inflation and black hole 
physics. The peculiar properties of the additional fields inducing dissipation is revealed by 
the covariantization. The links with the phenomenological approach to Quantum Gravity 
and with some Brane World scenarios are also discussed. 



^ This work has been presented in the SISSA conference From Quantum to Emergent Gravity in June 2007. 
It contains large extracts of the manuscript Confronting the trans-Planckian question of inflationary cosmology 
with dissipative effects. 
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1 Introduction 



Even though relativistic QFT provides an excellent description of particle physics, being non- 
compact, Lorentz symmetry has been tested only up to a certain energy scale pQ. Thus one 
cannot exclude the possibility that some high energy processes break the invariance under 
boosts, thereby introducing a threshold energy A^y, and a preferred frame. It is therefore of 
interest to understand what would be the signatures when this possibility is realized. 

Unlike dispersive effects which have been studied in details 12] , dissipative effects have so far 
received much less attention. In this work we provide the grounds for this extension. To this end, 
we first construct QFT displaying dissipative effects in the UV. Indeed, to handle dissipation 
requires settings which are wider than those of relativistic and dispersive theories because, if 
one introduces dissipation from the outset in the usual settings, one looses both unitarity and 
predictability. To preserve them, we shall therefore work with Hamiltonian theories in which 
dissipative effects are caused by interactions with additional degrees of freedom. Doing so, we 
shall discover that dissipative effects are generic. That is, when starting with a bare Lagrangian 
in which LI is broken by some kinetical or interaction term in the UV, the effective theory (the 
generating functional) unavoidably develops dissipation above a certain energy scale, simply 
because nothing can prevent this. (With relativistic QFT instead, LI did prevent it). We 
also learn that it is illegitimate to deal with purely dispersive QFT, since these do not satisfy 
Kramer's relations. 

Being engaged in a procedure of generalization, we should have a clear idea of our motiva- 
tions, aims, and requirements. We first present these aspects. 



1.1 Two requirements 

P*. Unitarity 

As already mentioned, we require that our QFT evolve unitarily. This requirement implies 
that the dissipative effects be produced by the Hamiltonian dynamics of the entire system. In 
other words, we shall introduce additional degrees of freedom called which play the role of 
an environment, and couple them to the original field (p in such a way that the latter develops 
dissipative effects. This means that the (dressed) two-point function of (p will be given by the 
usual QM trace 

Gwix,y) = Ti pT(i)ix)(j){y)], (1) 

where pt is the initial matrix density of the entire system + where 0(x) is the Heisenberg 
field operator evolved with the time ordered exponential of the total Hamiltonian, and where 
the trace is taken over both \1/ and 0. 

2'^. Stationarity and homogeneity 
Our second requirement concerns the properties of dissipative effects. When considering the 
theory in vacuum and in Minkowski space-time, we impose that the dissipative effects become 
significant only above a critical energy Alv, and that they preserve the stationarity, the homo- 
geneity, and isotropy of flat space-time. In this case, dissipative effects define a preferred frame 
which is inertial and globally defined. 
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Then, irrespectively of the properties of the additional degrees of freedom \& and their 
interactions with 0, the Fourier transform of the retarded Green function Gji{x,y), which is 
given by Gr{x, y) = 9{tx — ty) 2ImG'iy(x, y), where Gw is the Wightman function of eq. ([H), is 
of the form 

Gr{uj.p) = 7 ^ . (2) 

In the true vacuum, at the level of the 2-pt functions, the dissipative (dispersive) effects are 
indeed completely characterized by the imaginary and odd (real and even) part in uj of the 
(retarded) self-energy Sij(a;,p). 

In these expressions, the energy uj and the spatial momentum square p^ have been defined in 
the preferred frame. To prepare the covariantization of our theory, and therefore its extension 
to curved space-time, it is usefull to characterize the preferred frame in a coordinate invariant 
way by a unit time-like vector field, here after called l^. Then uo and p^ are given by 

uj = l^p^, p'=l.^'p^p,, (3) 

where 1.^^'^= r]^^ + l^l'^ is the (positive definite) metric in the spatial sections orthogonal to l^. 

The novelty is that is a function of uj and p separately, and not only of the relativistic 
invariant uj"^ — p^ as it is the case in relativistic QFT. When E depends on both uj and p, 
dissipation can become significant above a critical energy on the mass shell, i.e. along the 
minima of the denominator of eq. ([2]); a possibility forbidden in LI theories. 

With the observation of the simplicity of eq. ([2]), we understand that most of the properties 
of \E' will be irrelevant when restricting attention to observables built with only 0. In other 
words, the effective action of only retains little information about \E' in its moments, the second 
of which being the self-energy E. In this paper we shall exploit this freedom by choosing the 
simplest models of \E' which deliver the required properties of S. 

1.2 Aims and Motivations 

1.2.1 The phenomenology and the links with Quantum Gravity 

We first aim to describe the phenomenology of dissipative effects in Minkowski space-time and 
in vacuo. This is rather easy. From eq. ([2]) we understand that the phenomenology of dissipative 
effects respecting our (minimal) restrictions will be governed by a limited set of functions. In 
fact, at the level of 2-pt functions, the imaginary part of S is the only relevant quantity. One 
easily verifies that the following self-energies induce significant dissipation only above A^y 

Later in the text, we shall provide Lagrangians of \I' and (p which produce this class of imaginary 
self-energies. We can already relate this class to the set of non-linear dispersion relations which 
have been used in phenomenological studies: 

(2 \ 
^) +o{ipr'-')- (5) 
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The + (— ) sign gives super luminous (subluminous) propagation. It is also worth noticing that 
a phenomenological study of dispersive and dissipative effects taken toghether can be done with 
eq. (E]) by considering both ReS/j and ImSij non zero. 

Even though Quantum Gravity is not part of our main concerns, let us say a few words about 
the relations with our work. 

First, one should clearly distinguish two possible scenarios in which LI might be broken. 
First, the violation could be fundamental. In this case, the violation will show up in every cir- 
cumstances, including in Minkowski space-time (or what replaces it) and in vacuum. Secondly, 
the violation could be only effective. That is, it could be induced by radiative corrections in 
some circumstances, e.g. in black hole geometries, in cosmological backgrounds, or in thermal 
states. This line of thought have been developped in [31 S] 

Secondly, Quantum Gravity, whatever version is adopted, implies that the smooth manifold 
of GR should be replaced by a new structure when reaching a UV threshold energy, call it 
A^y. When adopting a phenomenological point of view, there is a shift in the interest. The 
question is no longer: What is this new structure ? but rather: How would it manifest itself 
in observables ? That is to say: What are the new expressions of the n-point correlation 
functions ? Remember that the predictions of a QFT are all based on its correlation functions. 

The simplest of these are the 2-point functions. In whatever replaces the Minkowski mani- 
fold, spatial homogeneity and stationarity will be preserved in the mean, (otherwise that version 
of QG would not describe our world). This implies that the Fourier transform of the "true" 
Feynman Green functioiJl will behave as in eq. ©I (Of course LI might still be respected in 
which case the self-energy induced by QG will only be a function of uj"^ — P^-) 

The link between our approach and the phenomenological approach to QG which consists 
in parametrizing its effects (rather than computing them from first principles) is clear: Since 
we provide the general expression of the 2pt function compatible with QM, in the sense that 

^Whatever QG may be, this function will be given by an expression similar to eq. ([1]) with the trace taken over 
the "true" degrees of freedom. For the benefit of ttie reader, I add fiere a remark made by a referee concerning 
ttiis sentence and my reply to it: "While of course one might claim this is a plausible scenario, one might as 
well conjecture that in no way transplanckian physics can be described by a quantum field theory of any sort 
(e.g. QM could be emergent as well)." My answer is twofold. First, I wrote "similar" precisely to leave open the 
possibility that the trace be not "the usual QM" one. Second, even if QM is only emergent, when considering the 
2-pt function of accessible (effective) degrees of freedom, as 4> in the present case, (or like the 2-pt function of 
the phonon field in a BEC), the true expression will, at least, contain a trace (to average out the configurations 
so that the correlation function only depends on the arguments of the two operators), a state function (to weigh 
these configurations), and two operators (which might be complicated composite objects), because, in the IR, the 
2-pt function must behave according to QM. The beauty of 2-pt correlation functions is that they are c-number 
functions depending only uj and p, both in the effective description and in the true description. Therefore a 
comparison of two versions directly delivers the relevant modifications, see for a comparison of 2pt functions 
obtained from QFT in a fixed background and from solutions of the Wheeler- DeWitt equation. 

^1% has been claimed that the reparametrization invariance of QG would severly restrict the true observables 
(e.g. to be topological charges). If this could be correct in pure 2+1 QG[6j, it is harder to conceive the 
relevance of this statement when applied to 4D, both from a phenomenological and a theoretical point of view. 
At a theoretical level, the quantum transitions of (heavy) atoms exchanging photons are governed by the 2-pt 
functions of the radiation field evaluated where the atoms sit, given their wave function. This should remain 
true when including QG effects, see [7j for a description of atomic transitions based on solutions of the WDW 
equation, i.e. when reparametrization invariance has been taken into account. 
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the Equal Time Commutations relations are satisfied, our expressions can also be used in a 
phenomenological QG perspective. 

1.2.2 Mode creation in expanding universes 

When assuming that LI is broken in the UV in Minkowski space, one can then also assume that 
the density of degrees of freedom, or modes, is finite and of the order of A|y. (This density 
should not be confused with the density of quanta which is always finite in decent quantum 
states, the "Hadamard" statesjU [9].) This option was not available as long as LI prevailed, 
and was in fact the source of the UV divergences. If the cutoff Klv is much higher than the 
typical frequencies involved in the observables we have access to, it is rather easy to show that 
the expectation value of these observables will not be significantly affected by A^y, a simple 
example which can be worked out explicitely is furnished by the Casimir effect. Therefore, 
in Minkowski space, one cannot expect any significant deviation induced by having cutoff the 
mode density. 

On the contrary, it is more challenging to consider how would such a theory behave in 
an expanding universe. Indeed one faces an alternative. Either the total number of degrees 
of freedom in a given comoving volume would stay fixed (or nearly fixed), and therefore the 
density would decrease like the inverse of the proper volume, or it is the density that would 
stay (nearly) fixed, and in this case, the number of degrees of freedom would grow linearly with 
the proper volume, thereby implying mode creation. 

The first alternative seems already excluded because the volume of our visible universe 
increased at least by a factor of e^^° since the onset of inflation (360 = 60 e-foldings during 
inflation x two for the radiation era x three because of 3D). Indeed either the density was 
absurdly high at the onset of inflation, or if the intial density was decent, we should today be 
lacking degrees of freedom. Moreover if it were true, one could measure the growth of the scale 
factor by probing locally the vacuum, thereby violating the Equivalence Principle. So, once 
having assumed a finite density in flat space, we are left with the conclusion that mode creation 
is unavoidable in expanding universes, and therefore in any curved background geometry. 

When adopting this second alternative, two questions should be confronted: 
How to describe mode creation in QFT ? 
What fixes the state of the newly born modes ? 

Our motivation is to confront these two questions in the presence of dissipative effects. To 
this end we first need to extend our QFTs to curved space times. 

The main principle we adopt is the Equivalence Principle (or more precisely its extension in 
the presence of the vector field l^). That is, the Lagrangians will be a sum of scalar functions 
of the four local fields 0, \E', g^y and which reduce locally to their Minkowski value in the 
zero curvature limit. This principle fixes the action density (up to the possibility of some non- 
minimal coupling) and determines the tensorial nature of the ^ fields. In the models we shall 
use, the ^ fields form a dense set of local degrees of freedom at rest with . 

When the Equivalence Principle is respected, without fine tuning nor additional hypothesis, 
we shall see that dissipative QFT are such that, as the universe expands, underdamped modes 
emerge from overdamped modes -thereby effectively describing the creation of propagating 
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degrees of freedom- in such a way that their (proper) density stays constant. Moreover, unhke 
in former attempts to describe "mode creation" [101 CI!; there is no need to supply an extra 
condition to fix the state of the newly born modes. Indeed, since our models are Hamiltonian, 
all the information is contained in the initial density matrix pT- It should also be noticed that 
when considering contracting universes, or locally contracting regions, unitary models imply 
the inverse process, i.e. "mode destruction", whereas it is unclear what the settings of [TT] 
predict in this case. As a last comment, the mechanism at work in unitary models can also be 
described as "mode conversion" by analogy to processes encountered in condensed matter]^ 

1.2.3 The trans-PIanckian question 

In inflationary cosmology, the primordial density fluctuations arise from the amplification of 
vacuum fluctuations which had very short wave lengths (very large proper frequencies) at the 
onset of inflation [12] . Similarly, Hawking radiation emitted by a black hole emerges from 
configurations which had extremely high initial frequencies [I3]. In both cases, the unbounded 
frequency growth questions the validity of the predictions because these have been obtained 
using the standard treatment, namely QFT in a curved space time. However there is no reason 
to believe that these settings still provide a reliable approximation for frequencies way above 
the Planck scale. 

Following original work of Unruh and Jacobsonp^ [15], the robustness of the standard 
predictions against modifying the theory in the UV have been tested by introducing dispersion 
relations which become non linear above a certain UV scale A^y, see eq. (|5]). Even though the 
propagation of the configurations is severely modified when this scale is reached, it was shown 
that the properties are essentially unmodified when the two relevant scales are well separated, 
i.e. in BH physics [T6] when k/Klv 1 where /t is the surface gravity, and in inflation [T71 ITS] 
when H/Klv <C 1 where H is the Hubble parameter. What guarantees the robustness is that 
the vacuum state evolves adiabatically. 

Our aim is to generalize these works by providing dissipative models in which the power 
spectrum (and Hawking radiation) can be computed. These models are the same as those of 
the former subsection. The fact that the Equivalence Principle is preserved will guarantee, as 
we shall see, the adiabaticity of the evolution of the true vacuum as long as the gradients of 
the metric are much smaller than the UV scale 

To conclude, we would like bring to the reader's attention to the following remark. To 
obtain the power spectrum in inflationary models, and the asymptotic properties of Hawking 

^To obtain dissipative QFT, we could have searched for inspiration in quantum optics or other condensed 
matter models since the effective LI is also broken above a certain scale. We have chosen not to pursue this 
approach for several reasons. First this approach lacks generality and thus hides the steps necessary to construct 
dissipative models. The second reason comes from the Equivalence Principle considered in expanding universes. 
Indeed any QFT emerging from a discrete structure, e.g. a set of atoms, is bound to violate the Equivalence 
Principle after a large number of e-folds. We see no way to escape this conclusion. 
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radiation, it is always sufficient to know the 2pt function of eq. ([T]) since 




dte'"^^ Ga(t,r;0,r,l,m), r Th. 



'xe~'''''Ga{t,x;t, 0), 



(6) 



In inflation, the first equality follows from the definition of the power spectrum Pk which is 
given by the spatial Fourier transform of the anti-commutator Ga evaluated, at equal time after 
horizon exit, in the Bunch-Davies vacuum. As of Hawking radiation, n{u!, I, m), the asymptotic 
distribution at fixed angular momentum l,m, is given by Fourier transform with respect to the 
asymptotic time of the in-vacuum anti-commutator far away from the hole. 

We now remark that the knowledge of Ga becomes necessary in the presence of dissipation 
because the only way to obtain Pk and n{uj) is through the above expressions. 

We have organized the paper as follows. We first construct dissipative theories in Minkowski 
space-time. In this case, because of stationarity the analysis is simple and instructive. We then 
covariantize these models and briefly comment on the trans-Planckian question in inflation 
and in black hole physics. In long Appendices we provide self-contained presentations of the 
properties of 2-pt functions in the presence of dissipation. 

2 Dissipation in Minkowski space from LV effects 

In this section, we provide a class of models defined in Minkowski space time which exhibit 
dissipative effects above a certain energy scale A^v Stationarity, homogeneity and isotropy 
will be exactly preserved. Therefore, the only invariance of relativistc QFT which is broken is 
that under boosts. These theories define a preferred rest frame which is globally defined, as it 
is the case of FLRW space-times. Even though a covariant description exists, for simplicity of 
the presentation, we first work in that frame. At the end of this Section we shall covariantize 
them to prepare the extension to curved space-times. 

2.1 Free field settings 

We start with a brief presentation of the free field description to introduce our notations and 
to point out what are the properties which are lost in the presence of dissipation. 
The action of our free massless field (p is the usual one: 



Due to the homogeneity of space, the equation of motion can be analyzed mode by mode: 




(7) 




(8) 
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The mode operator (f)p{t), the Fourier transform of the field operator, obeys 

{d^ + coD^p = , (9) 

where cj^ = = p ■ p is the standard relativistic dispersion relation. Notice that eq. is 
second order, homogeneous (no source term), and time reversible (no odd power of dt), three 
properties we shall loose when introducing interactions breaking LI. 

In homogeneous space-times, the canonical Equal Time Commutator of the field and its 
momentum implies 

[0p(t),a4,(t)]=z5=^(p-p'). (10) 
When decomposing the mode operator as 

0p(t) = «p0p(t) + alp0;(t), (11) 

where the destruction and creation operators satisfy the usual commutators 

[ap,aj,,] = 5^(p - p') , [ap,ap/]=0, (12) 

eq. fllUp is verified at all times because the Wronskian of the positive frequency (c-number) 
mode 

Mt)=e~'^^'/i2u;,y/\ (13) 

is constant (and conventionally taken to be unity). 

Had an odd term like 'ydt been present in eq. (Q the constancy of the Wronskian would 
have been lost. Hence the possibility of realizing the ETC fllOj) with the help of eq. 0121) would 
have been lost as well. This already indicates that, unlike dispersive (real) effects, dissipative 
effects require more general settings than the above. 



2.2 Interacting models breaking LI, general properties 

We now introduce additional degrees of freedom, here after collectively named which induce 
dissipation above the energy A^y. We shall work with a particular class of simple models in 
order to get an exact (non-perturbative) expression for the two-point function of of eq. ([T]). 
Before introducing these models, we derive general results valid for a// unitary QFT's possessing 
dissipative effects above Alv in the ground state (the interacting vacuum). 
We assume that the total action decomposes as 

St = + + Sff,^^ , (14) 

where the first action is that of eq. ([7]), the second one governs the evolution of the \1/ fields, 
and the last one the coupling between and these new fields. 

We also impose that the last two actions (and the state of the system, the density matrix 
Pt) preserve the homogeneity and isotropy of Minkowski space. From now on, the cartesian 
coordinates t, x are at rest with respect to the preferred frame which is defined by the action 
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5*^ + 5*^,*, i.e. in a covariant notation dt = I'^d^, see eq. (E]). For these models, the Fourier 
transform of the Wightman function of eq. ([T]) is of the form 



Gp,p'(t,t') 



Tr[pT0p(t) 0^,(0], 



(15) 



At this point, an important remark should be made. In the presence of interactions, the notion 
(and the usefulness) of the time-dependent modes given in eq. (fT^ disappears. Instead, the 
time-dependent function Gwit,t';p) of eq. (ITSl) is always well-defined, for "all" choices of \E' 
and for all actions S^j + S*,^,*. We shall return to this point in Appendix B after eq. (1791) . 

When the situation is stationary, Gw{p',t,t') further simplifies in the frequency representa- 
tion: 



When working in the ground state, one reaches the simplest case in the following sense. What- 
ever the fields may be, the Fourier transform of the time ordered (Feynman) propagator, is 
always of the form given in eq. ([2]), and therefore characterized by a single function, the Feyn- 
man self-energy T,p{uj,p). When restricting attention to Gaussian models, T,p{u!,p) is given 
by a 1-loop calculation, whereas it contains a series of IPI graphs for non-Gaussian models. 
In non-vacuum states and in non-stationary situations, 2pt functions and self-energies have a 
more complicated structure p?j. To understand this strucure it is useful to study separately 
the commutator Gc of 0p, which is odd in u, and the anti-commutator Ga, which is even. 

Let us conclude with two remarks. First, the effective dispersion relation of is a posteriori 
defined by the poles of eq. ([2]). In this way, non-trivial dispersion relations arise from dynamical 
processes rather than from being introduced from the outset. The present work therefore 
provides physical foundations (and restrictions, as later discussed) to the kinematical approach 
which has been adopted in the literature. Second, from analyzing dynamical models, we shall 
see that, even in the vacuum, on-shell dissipative effects (i.e. dissipation arising along the 
minima of the denominator of eq. ([2])) are unavoidable when LI is broken in the UV by the 
action S*^ + S*,^,*, in complete opposition with the fact that on-shell dissipation is forbidden 
when working (in the vacuum) with LI actions. 

2.3 Gaussian models 

To simplify the calculation of S(ti;,p) and to get non perturbative expressions, we assume that 
the action St is quadratic in all field variables. At first sight this could be considered as an 
unjustified hypothesis. However, it should be recalled that we are not after computing S from 
first principles. Rather we aim to compute the power spectra of ([6]) given the properties of S, 
and this can be done with quadratic models. 

Since we are preserving the homogeneity of Minkowski space, the Gaussian assumption 
implies that the total action splits as 




(16) 




(17) 
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where each action St{p) depends only on the complex mode operators </)p, ^I/i(p). (Since 
(f)^p = (p-p, ^i(p) = p), St{p) + St{—p) is real.) The structure of St{p) is 



where i is a discrete (or continuous) index, where fli{p) is the frequency of the oscillators ^'i(p), 
and where gi{p) is the coupling constant at fixed p,iu 

Models of this type have been used since many years (at least since the early 50's) and 
for different purposes. They have been introduced (in their continuous version) to study non- 
pertubatively atomic transitions and to model quantum electrodynamics in the dipole approxi- 
mation, see [22] and refs. therein, see also [23l [2l] for an application to the Unruh effect. They 
have been used in Quantum Optics [25], to model quantum Brownian motion, and to study 
decoherence effects [23. Needless to say that our intention is not to provide a new approach 



Instead of chosing a priori a specific model, we shall solve the equations of motion without 
specifying the set of \I/i(p), their energy Qi{p) and their coupling gi{p). We shall choose them 
in the next subsection to further simplify the equations we shall need to solve. To preserve 
stationarity in Minkowski space, the fi's and the g^s must be time independent. When ^f{p) 7^ 
Mf +p'^, the kinetic action of ^I/i(p) breaks LI and defines the preferred frame. On the contrary 
when fif {p) = Mf + p^ the preferred frame is only defined by the interaction term through 
the p-dependence of the coupling functions gi{p). 

The equations of motion are 



^When this action was presented in seminars, it was often found ad hoc and/or artificial. We dispute both 
charges for the following reasons. Concerning the ad hocery, it should be realized that any unitary (homogeneous) 
model exhibiting dissipation will be governed by an action of the type pH]). either strictly speaking, or in the 
sense that the 2pt function of eq. ([T]), will exactly behave as that obtained from a model governed by an action 
of the type (|18|) . Given this, the artificial qualification appears to us also unjustified since the 2-pt functions 
derived from "physically well-founded" models would perforce belong to the class derived from eq. pS)) . It 
should be noticed that the ad hoc qualification has been also often used against dispersive models, apparantly 
not realizing that the robusteness of the spectra was established essentially for all dispersion relations [TSl [Mj . 
What was probably meant by is that it is not yet known which of these relations is actually revelant, but this 
is another issue. 

^Depending on the point of view adopted, they can be solved and analyzed by means of different methods. 
In what follows we shall use the simplest approach based on Heisenberg equations of motion. Since we shall 
work at zero temperature, this approach is appropriate because quantum aspects dominate over stochastic ones. 
Even though legitimate, we shall not use the general methods [25] which have been developed to study "open 
quantum systems" since these methods (Influence Functional, Master Equation, Closed Time Path Integral) 
somehow hide the simplicity of the models we are dealing with. In fact only but standard Quantum Mechanics 
is required to solve them. Moreover we are planning (in a subsequent work) to study the correlations between 
(j> and '^i. Therefore we shall treat cj) and '^i on equal footing as in [5^ 1251 [51] . 





(19) 
(20) 
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The general solution of the second equation reads 

v[/,(p,t) = vl/°(p,t) + j dt'R1{t,t'-p)g,{t';p)Mt'), (21) 

where \E'°(p, t) is a free solution which depends on initial conditions imposed on \E'j(p). The 
second term contains R°{t,t';p), the (free) retarded Green function of ^E'j(p). It obeys 

{d^ + QUp))RKt,t';p) = 6{t-t'), (22) 

and vanishes for t < t'. (To prepare the application to time dependent geometries, we have 
treated Up, flf and gi as arbitrary functions of time (in cosmology these become indeed time 
dependent through their dependence in the scale factor a(t)). ) Injecting eq. ( 12T1) in eq. ( |T9i) 
one gets 

{d^ + ijl)(j)p = E,g,{t;p)^1{p,t) + E,g,{t;p) J dt' R°{t,t';p)g,{t';p)^p{t') . (23) 

The (exact) solution of this equation has always the following structure 

Mt) = + / dt'Gr{t,t'-p)[J:,g,{t';p)^''^{p,t')] . (24) 

The first term is the "decaying" solution. It contains all the information about the initial 
condition of (pp, and obeys the non local equation 

j dti [ 6it - ti){dl +ujl)-T.i g,{t- p) R1{t, ti-p) (7.(ti; p)]0p(ti) = . (25) 

The second term is the "driven" solution. It is governed by the initial conditions of ^E'°(p) 
and by the (dressed) retarded Green function, the solution of eq. (1^ with 5{t — ti) on the 
r.h.s. Therefore the evolution of both 0*^ and Gr fully takes into account, through the non-local 
term in the above bracket, the backreaction due to the coupling to the additional degrees of 
freedom. In Gaussian models, the backreaction is quadratic in gi. Hence the solutions (f)'^ and 
Or are series containing all even powers of gi. Therefore Gaussian models do give rise to non 
perturbative effects. Moreover since gi{t;p) are arbitrary functions of p and t, at this point 
there is no reason to consider non-Gaussian models. 

We conclude this subsection with two remarks. First, eq. ([21]) also furnishes the exact 
solution for the Heisenberg operator 0p because the equations we solved are all linear. Since 
we work quantum mechanically, it is necessary to study the 2pt functions of 0p. We review 
their properties in the Appendices. 

Second, eq. (P^jl tells us that in general, the coupling to an environment (even a linear one) 
gives rise to non-local equations of motion. When dealing with stationary situations this does 
not cause any problem because, as shown in Appendix B, observables can be (algebraically) 
computed in the frequency representation. However in non stationary situations and in curved 
space-times, to be able to compute observables it becomes imperative to simplify eq. (l25i) . A 
question we now address. 
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2.4 Time dependent settings 

In Appendix C, we provide a class of stationary models characterized by the power of the 
ratio p/A^y which specifies how dissipative effects grow with the energy, see eqs. (jH IHTj) . 
This class covers the general case and can be used as a template to study the phenomonological 
consequences of dissipative effects. In addition, after eq. fl7^ . it was noticed that the dissipative 
properties are governed by the product g'^{uj,p)R°{u!,p), as can be understood from eq. ( l25l) . 
Therefore theories with different degrees of freedom \l/ and with different coupling terms but 
delivering the same product will give rise to the same (stationary) phenomenology. 

In this Section we exploit this freedom to simplify the expressions having in mind the 
transposition of our model from Minkowski space to curved metrics. Therefore the models we 
shall use should possess two-point functions with simple properties when expressed in the terms 
of space-time coordinates (and not in Fourier components uj,p). 

The core of the problem is that, when dissipation is strong, the relevant observable Ga, 
see eq. is given by the "driven" term in eq. (lU^ . which is a double integral containing 
a kernel (which can be computed) and the retarded Green function of which is not known. 
Indeed it is only implicitely defined as a solution of eq. f l25l) . which is, in general, non-local. 
We are thus led to choose the action of \E' in order for eq. ( l25i) to be local. Notice that this 
requirement concerns the Green function of the \Ef-environment and not that of the 0-system 
we probe. Thus it will not restrict the latter's phenomenology. 

Given this choice, an appropriate class of models is defined by the action 



+1 jdt jdk ^*{p,k){-d^ - {nALykf)^{p,k) 

n+1 



+gALv Jdt Jdk (^-^^ 0p dt^*{p, k) . (26) 

When compared with the action of eq. (ITS!) . we have replaced the discrete index i by the 
integral over the dimensionless variable k. As recalled in Appendix B, one needs a continuous 
spectrum of the environment to have proper dissipation, see discussion after eq. (1721) . The 
continuous variable k can be viewed as a momentum (expressed in units of A/^y) in an extra 
fiat fifth spatial dimension. The relationship with the Brane World Scenarios of [20] is clear. 
In the 'atomic' version of this model [231 121] which has inspired us, the radiation field is a 
massless 2 dimensional field propagating in t and in the dimension associated with k. 

In Sfi)i^ we have factorized out A^y so that the coupling constant g is dimensionless. We 
have also introduced an additional time derivative acting on \E'. This choice leads to the above 
mentioned 5{t — t'). Indeed, on the one hand, taking this extra derivative, the continuous 
character of k, and the fact that g is independent of k, eq. ( I23l) becomes 

id^+ujl)(l)p = gndt Jdk^\p,k,t)- g^dt J dt' JdkR\t,t';k,p)dt\gn<P^{t')), (27) 

where gn = gA.Lv {p/-^lv)"'~^^ ■ On the other hand, for each 3-momentum p, the Green function 
of S = J dk'^{k) is that of a massless 2-dimensional free field. In Fourier it is given by 
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R°{uj, k) = l/{-{uJ + ief + {nKLvkf)- Hence it obeys 



dtW{t,t') = dt 



I 




dkR"{uj,k) e 



6{t - 1') 



(28) 



which is the required property to simphfy eq. (!2711 . 

When is constant, the retarded Green function of associated to eq. (|27|) obeys the 
following local equation 



To make contact with Appendix B and C, let us rewrite this equation in Fourier transform, 



We thus see that ReS,. = and that ImS^ is (exactly) given by times the expression of 
eq. ( ISTj) . Thus, even though we have chosen a simple and definite form for R°(t, t'), the above 
action does deliver the n dissipative behaviors of Appendix B by choosing the appropriate 
power of p in the coupling gn appearing in the action S(f,if,. 

When g and are arbitrary time-dependent functions, the Fourier analysis looses its power. 
However, in time dependent settings, our Gr still obeys a local equation: 



where the n-th decay rate 7n(i) = g'^{t)ln is now a definite time dependent function. When 
transposing the model of eq. fl26|) in expanding homogeneous universes, the Green function will 
therefore also obey a local equation. 

2.5 Covariant description 

We now provide the covariant version of the action of eq. This expression will then be 

used to define our theory in curved backgroundslll 

^Perhaps the motivations for these two steps require further explanation since another referee wrote: "/ do not 
see the point of a covariant description of the \1/ field since it is nothing more than a convenient parametrization 
of the environment degrees of freedom." Several points should be mentionned. First, as in studies of dispersive 
effects, we require that our QFT be well defined in arbitrary backgrounds. It would be physically questionable 
(or even meaningless) to study dissipative effects on the observables of eq. © from theories which can be defined 
only in a particular class of space-times. Second, we have introduced the ^' field to parametrize the effects of 
dissipation which would stem not from any theory, but only from fundamental theories obeying the Equivalence 
Principle, see jTS] for a prototype. Third, from the point of view of open quantum systems, it could indeed seem 
inappropriate to proceed to a covariantization, since in most situations (heat bath, condensed matter systems) 
there is a preferred frame which is globally defined. However, when the system is non-homogeneous (e.g. a fluid 
characterized by a non- homogeneous flow), low energy fluctuations effectively live in a curved geometry [29j. 
Moreover, in this case, short distance effects, i.e. dispersive (or dissipative) effects, are covariantly described 
when using this metric [TH [H] , because they arise locally. In brief, we shall write down actions in curved space 
using the EP which requires a covariant description. 




(29) 




(30) 




(31) 
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Two steps should be done. We need to go from p considerations to a local description, 
and express the various actions in terms of the unit vector field F and the spatial metric -L^'^ 
defined in eq. ((31). Both are straightforward and, in arbitrary coordinates, the total action 



S, 



(n) 



Jd^pS, 



T ' 



S. 



in) 



p) reads (in Minkowski space-time) 
1 



d Xy 



1 

+ 2 



d^X\/—g / dk 



LV 



d X\ 



-9 



A 



Cv|> 



{n+l)/2 



Fd„ dk^ 



(32) 



where A is the Laplacian on the three surfaces orthogonal to F. 

We have slightly generalized the action of eq. (1261) by subtracting c|, -L'^'^ to the kinetic term 
of \E', where c|, <C 1. With this new term, p-th components of the field are now massive fields 
which propagate with a velocity (relative to F) whose square is bounded by c|,. In addition 
they now possess a well defined energy-momentum tensor which can be obtained by varying 
their action with respect to g^'^. To obtain the simplified expressions we use in the paper the 
(regular) limit c|, — should be taken. 

In this limit, the (free) retarded Green function of the S = J dk'^ field obeys a particularly 
simple equation when expressed in space-time coordinates. Indeed, it obeys 

1 5^(x'^-?/^) 



A 



LV 



-9 



(33) 



On the r.h.s, one finds the "invariant" delta function with respect to the invariant measure 
d^x^—g. This equation is nothing but the covariantized and "localized" version of eq. ( l28l) . 
Its physical meaning is clear: the back-reaction on 0(x) through S is strictly local. 

Nevertheless, the equations of motions do not have a particularly simple form in space-time 
coordinates. Using the condensed notation S = J dk'^ ^ the relevant equations read 



.d,4^gg^'d,. 



\x) 



\/-9 

where the interacting field S field is 



9^ 



LV- 



9 ^^LV 



n+1 
2 



gALv / d^yV^ R°(a;, y) 



:d„ I 



A 



, A2 

^^LV 



n + l 



(34) 



(35) 



When inserting eq. (l35l) in eq. ( IMl) . using eq. ( l33l) . one verifies that the dissipative term is 
local and first order in Fdf^. Dissipation thus occurs along the preferred direction. 



3 Dissipative effects in curved space-time 
3.1 The Lagrangian and the properties of ^ 

To define a dissipative QFT in an arbitrary curved geometry, we need some principles. From a 
physical point of view, we adopt the Equivalence Principle, or better what can be considered 
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as its generalization in the presence of the unit time-hke vector field. We are in fact dealing 
with two (sets of) dynamical fields, the field we probe, and the fields we do not; but also 
with two background fields (7^1, and l^. The Generalized Equivalence Principle means that the 
action densities of the dynamical fields be given by scalar functions (under general coordinate 
transformations) which reduce to those one had in Minkowski space time and for a homogeneous 
and static l'^ field, i.e. those of eq. ( l32l) . 

The densities are thus not completely fixed by the GEP, as it was the case with the EP. 
For a scalar field in a curved geometry, there was always the possibility of considering a non- 
minimal coupling to gravity by adding to the lagrangian a term proportional to In the 
present case, the ambiguity is larger because the vector l'^ allows to form new scalars, such as 
its expansion 

&^V,F = {-gr'^%[{-gy/'n, (36) 

where is the covariant derivative with respect to the metric gfj,^. The ambiguity can only 
be resolved by adopting some additional principle, such as the principle of minimal couplings 
which forbids adding such scalars. 

Rather than adopting it, our second principle (which simplifies the equations of motion, but 
which is by no means necessary) is that the locality of the back-reaction effects of (p through 
\1/ be preserved. That is, we choose the non-minimal coupling of ^ so as to keep eq. ( I33l) . 
Starting from eq. ( l32l) . the locality is preserved by replacing both in S"^ and Sifj^f, 

Pa^^, ^ Vi^Hk = I'^d^'^k + I = ^ (/'^V^^fe + r^k]) ■ (37) 

To simplify the forthcoming equations, we use the fact that one can always work in "pre- 
ferred" coordinate systems in which the shift vanishes and in which the preferred time is 
such that Iq = 1. (We assume that the set of orbits of l'^ is complete and without caustic. In 
this case, every point of the manifold is reached by one orbit.) It is useful to work in these 
coordinate systems with rescaled fields = {—gY^^"^ because one can then group the above 
two terms into a single expression: 

Vi^k = i-g)-'/' I'd.ii-gfl^ VP,) = {-g)''/' Fd,^>l . (38) 

Notice also that there exists a subclass of background fields {g-,1)-, for which one can find 
coordinate systems such that both the shift and g°^ vanish. In these comoving coordinate 
systems, the above equations further simplify since only the spatial part of the metric matters 
because —g = he where h = det(±ij). 

Having chosen this non-minimal coupling, one verifies that the kinetic term of the rescaled 
fields is insensitive to the "curvature" of both (7^,^ and l'^ (when the limit c|, — > is taken). 

^ To illustrate this point, let us consider FLRW flat metrics. In comoving coordinates, one has ds^ — 
—dt^ + a^dx^, and h = he — a^. In Lemaitre coordinates X — ax, one has ds^ = —dt'^ + {dX — Vdt)^, where 
the velocity is V = HX. The spatial sections are now the Euclidean space with h ~ 1, and the (contravariant) 
components of the unit vector field are = 1, V. To compute he one should solve the equation of motion of 
comoving (free falling) observers dX — Vdt = 0, and use the initial position as new coordinates. This procedure 
will be done in subsection 13.31 starting with Painleve-GuUstrand coordinates to describe the black hole metric 
and P field. 
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Moreover the differential operator which acts on the retarded Green function of in the 
equation of motion of 0, see eqs. ( I34l 1351) . is also "flat" thereby guaranteeing that the modified 
version eq. ( |33l) still applies, that is 



where Ii°{x,y) is the retarded Green function of the rescaled field It obeys (in pre- 

ferred coordinate systems) dfH" = ^^^/A^y, and "defines" the retarded Green function R° = 
(— (7)^/^R° {—gY^^ which is a bi-scalar. Hence the equations of motion in an arbitrary back- 
ground " tensor- vector metric" specified by the couple {g^uj'^), are given by eqs. ( IMl 1351) with 
the substitution of eq. (!37|) . 

Several remarks should be made. First, from the simplified equation 9jR° = 6'^/Alv it 
might seem that the background tensor metric gfj^^, plays no role. This is not true, it enters 
indirectly as it is used to normalize the vector field at every point. Therefore, it intervenes in 
the specification that the (proper) energy scale A^y stays fixed as the universe (or the comoving 
volume he) expands (as required by the Equivalence Principle). 

Second, in the limit c|, — >■ 0, the (rescaled) fields define a new kind of field. They 
propagate in an effective space-time given by the time development of the 3-dimensional set 
of orbits of the field. Indeed, at fixed k, \l/fc(a;) can be decomposed in non-interacting local 
oscillators, each of them evolving separately along its orbit. This situation is similar to the 
long wave length (gradient-free) limit of cosmological perturbations considered by Steward and 
Salopek[3n]. In the absence of l^, the space time geometry must be (nearly) homogeneous for 
the mode action to possess this decomposition. However, when the l'^ field is given, one can 
identify, even in non-homogeneous metrics, every point in space-time in an invariant way by the 
spatial position of the corresponding /'^-orbit at some time, and the proper time along the orbit 
(as long as has no caustic). We can thus build covariant actions exploiting this possibility 
and consider fields composed of a dense set of local oscillators (i.e. one has 1 degree of freedom 
per point) at rest with respect to Z'^. The fields belong to this class of fields, here after 
called Chelva fields. 

Third, the analogy between the above formalism and the DeWitt [28] way to handle the 
short distance behavior of Green (Hadamard) functions seems worth deepening. 

Fourth, even though we have chosen to simplify the equation of motion of by imposing 
that the rescaled kernel obeys c?iR° = ^^^/A^y, the 2pt functions of in curved backgrounds are 
highly non trivial, and in particular the anti-commutator Ga- Indeed Ga is given by a double 
time integral, see eq. (l65l) . governed by the kernel N which cannot be approximated by a delta 
function in the vacuum, see Section 3 in [3T] . 

3.2 Application to cosmology 

To get the equation of motion we consider eq. ( l32i) (with the curved metric modifications 
introduced in the former subsection) in a FLRW metric ds^ = —dt^ + a?{t) d'x.d'^. To simplify 



1 5^(x^ 



^/^ra,(R°(x,y)) i-giy))-'^' 



(39) 
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the notations, we use the conformal time drj = dt/a and work with the rescaled fields 0^ = acf) 
and '^r{k) = a^/^\l/(fc). Notice that their power in a differs. Dropping these r indices, working 
in Fourier transform with respect to the (dimensionless) comoving coordinates x, the equation 
of motion of Heisenberg operator (pp is 

+ 27„ d, + {ooliT^) + a,7„)) 0p = 9n a,S°(p) , (40) 

The conformal frequency ujp{i]) = — d^a/a is that of a rescaled minimal coupled massless 
field. In this expression, as everywhere in this subsection, p is the conformal (dimensionless 
and constant) wave vector. The time dependent coupling coefficient is 

gn = g a'/' Alv {p/aALvT^' ■ (41) 

Its unusual a-dependence follows from the different rescaling of the two fields. Straightfor- 
ward algebra gives 7n(^), the decay rate in conformal time. When compared to the comoving 
frequency p to get the relative strength of dissipation, one gets 

P 2 pAlv 2 aAiv 2 Alv 

In the last equality we have re- introduced the proper momentum Pphys{,o) = p/a. With this 
equation we verify that, at any time in an expanding universe and for any mode 0p, the relative 
strength of dissipation is simply that obtained in Minkowski space at the corresponding energy 
scale, see eq. (IS^ . This is nothing but the expression of the Generalized Equivalence Principle, 
at the level of the dynamical equations. 

Notice however that the equation of motion in an expanding universe contains a frequency 
shift 

a,7. = -(2. + l)/(^)(^)^-\ (43) 

We have factorized out to get the relative value of the shift. It vanishes both when dissipation 
is negligible and when aH/p -C 1, i.e. when the expansion rate H is negligible with respect to 
the proper momentum. From this expression we can already conclude that it cannot play any 
role when the two relevant scales H and Alv are well separated, i.e. when 

cr = -^«l. (44) 

Indeed when the physical momentum is high and of the order of Alv, the relative frequency 
shift is proportional to a, and when the physical momentum is of the order of H (at horizon 
exit, see the Figure), it is proportional to cr^"+^. 
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Figure caption. We have represented in a log-log plot and by a dashed line the evolution of 
dn = Rh/o, the Hubble radius in comoving coordinates, as a function of a, both during infla- 
tion wherein dn oc 1/a and during the radiation era wherein dn oc a. We have represented by a 
thick line the trajectory of the cutoff length scale c/a = l/aA^y in the case where l/A^y <^ Rh 
during inflation. The dotted line corresponds to an intermediate fixed proper length A which 
obeys l/A/^y-C A-C-R^'^''. The vertical line represents a fixed comoving scale dp = 1/p. Below 
the cutoff length, all modes are overdamped. When a mode exits the cutoff length, it becomes 
underdamped and starts propagating. When it crosses the intermediate length \, it behaves as 
a free mode, and gets amplified only when exiting the Hubble radius. As explained in the text, 
adiabaticity, which is guaranteed by 1/A.lv Rh, guarantees in turn that, near \, modes are 
all born in the Bunch-Davies vacuum when the environment ^ is in its ground state. 

Let us now establish that the power spectrum is robust, i.e. in spite of strong dissipative 
effects, when eq. (1441) is satisfied, the standard expression is recovered when the \1/ fields are in 
their ground state. 

In free settings, when assuming that all 0p are in their ground state at the onset of inflation, 
i.e. working in the Bunch-Davies vacuum, the anti-commutator, see eq. (IMl) . evaluated at equal 
time rj is simply given by 



In the last equality we worked for simplicity with a constant Hubble parameter if, in which 
case a = —l/Hrj. Then the power spectrum of the physical (un-rescaled) field given by 



becomes constant after horizon exit. We have added a p subscript to H because in slow 
roll inflation, the relevant value of H{t) for the p-mode is that evaluated at horizon exit, i.e. 
Hp = H{tp), where tp is given hj p/a = H. The above equation shows that Pp acquires some 




(45) 



Hi (1 + {pvf) 



(46) 
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scale dependence through Hp. Similarly the deviations from this standard behavior stemming 
from some UV modification of the theory will also depend on p through Hp (and its derivatives). 
For an explicit example, we refer to [35j where the modifications of the spectrum stem from 
the fact that the adiabatic vacuum is imposed when p/ain is taken large but finite. 

In the presence of dissipation, the expression for Ga radically differs from the above. Indeed, 
when dissipative effects grow with the energy (as we suppose it is the case), one reaches the 
conclusion that in inflation the decaying solution of eq. (HOj) is completely erased (unless one 
fine-tunes the number of e-foldings so as to keep a residual amplitude). That is, the mode 
operator is entirely given by its driven term, the second term in eq. i^^. Hence the power 
spectrum is also purely driven and given by the second term of eq. 



Gi~{v,v,P) = Jdvi Jdv2Gr{v,Vi,p)Gr{v,V2,p)^{vi,V2,p), (47) 

where Or is the retarded Green function, solution of eq. ( I40l) with 6(1] — rji) as a source, and 
where the kernel N is the anti-commutator of Qnl^d^'^p, the source of eq. (HOl) . From eq. (H71) 
we learn that only the (initial) quantum state of the environment matters. In other words, 
because of the strong dissipation at early times, the power spectrum is independent of the 
initial state of 0, whatever it was. 

In spite of these differences, when eq. (jS]) is satisfied and when the environment is in its 
ground state, the predictions are unchanged, i.e. the power spectrum obtained from G^™*^"- 
coincides with that obtained with G{^'^'^ because the combined evolution oi cj) + consists in a 
(adiabatic) succession of stationary states ordered by the scale factor a. 

The proof goes in two steps. First scale separation, guarantees the adiabaticity of the 
evolution. That is, when 77 and rj' are close (in the sense that 1 — a{ri)/a{ri') -C 1), Gf"™^^ is 
well approximated by 



G' 



driven 



M,p)^GT''\n-n'-u;p{a),gp{a)) = J ^ e''^^^""') G.(a;; u;,(a), ^7,(a)) , (48) 



where Ga{uj; u!p{a) , gp{a)) is the Fourier component of the anti-commutator calculated with the 
Hamiltonian characterized by the constant values of the frequency ujp{a) and the coupling gp{a) 
both evaluated with a = a{rj). In Appendix B, the stationary value of these anti-commutators 
has been algebraically solved for all frequencies and all couplings, see eq. ( !77|) . Moreover, since 
by hypothesis, we are in the vacuum, eq. ( 1751) also applies with n{uj) = 0. In other words, the 
value of the anti-commutator follows from that of Gc- This is enough as it guarantees that when 
the mode (pp becomes free, i.e. much after A^^y-exit but before horizon-exit {H <^p/a <^ A^y), 
eq. fl80p applies, see the Figure. Thus, irrespectively of what the coupling with the environment 
was, adiabaticity implies 

G'r^-iv. V,p) - Gf-(r?, V,P) = , (49) 

2ujp[a) 

thereby guaranteeing that no modification of the power spectrum will be found. In brief adi- 
abaticity means that the evolution proceeds slow enough for not inducing any (non-adiabatic) 
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transition, which in the present case would correspond to creation of pair of (j) quanta with 
opposite wave vector p. 

The second part of the proof consists in providing an upper bound for the probabihty 
amphtude of obtaining a non-adiabatic transition. This amphtude is governed by the relative 
frequency change 

where Pphys = pIcl and 7phi/s = 7/*^ are the physical momentum and decay rate. Therefore, 
going backwards in time from the free regime in the underdamped regime, up to p^^ys > ^Iphys 
(i.e. Pphys < A^y), the non adiabatic parameter raises from zero but stays bounded by 

d^UO'^^^ /rt I. \ 4n+l 

-^<3(2n + l)a(f^^) < 3(2n + 1) a « 1 . (51) 

This guarantees that the amplitude for the system to jump out the ground state is bounded 
by a (up to an overall factor which plays no role). There is no need to study the stability 
of the ground state in the transitory regime from underdamped to the overdamped modes 
(for Pphys > A^y), because whatever transitions happened is suppressed by a factor e"-^'^'^* ~ 
exp(— l/cr(2r?, + l)) ^ 1. This completes the proof that scale separation guarantees adiabaticity. 



3.3 Dissipation in Black Hole metric 

In spite of the fact that the background metric is no longer homogeneous, the effects of dissi- 
pation on Hawking radiation can be studied along lines similar to the above. 

For simplicity we consider only spherically symmetric and stationary BH metrics. We also 
choose the unit vector field /'^ to be stationary and associated with Freely Falling observers 
which start at rest at infinity. In this case, the expressions for both the metric and /'^ simplify 
using the PG (Painleve-Gullstrand) coordinates, t, r, 6, 0. One has 

ds'^ = -dt^ + {dr - vdtf + r'^d^'^ , 
= (l,t;(r),0,0), 

g^Lv = -lf,lu+^f,u, with ±^^,= diag(0, l,r^r^sin^6') , (52) 

where v{r) < is the radial velocity of the FF observers, and t their proper time, not to be 
confused with the " Schwarzschild" (Killing) time. From the last equation we learn that the 
3-surfaces perpendicular to are simply the Euclidean space. Therefore h = det -i-ij= sin^ 9 
is independent of the expansion of f , as it is the case in cosmology in Lemaitre coordinates, see 
footnote [HI 

However PG coordinates are not comoving along the /'^ field since the shift is given by v. 
To determine the space-time dependence of the comoving volume element y/hc, we introduce 
the "comoving" coordinate tq = ro(t,r) defined by 

r dr' 

Jro v{r') ^' ^^^^ 
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By definition r = r{t,ro) gives the trajectory of tlie FF observer wfio started from ro at t = 0. 
Tfius, for every {t, r), tq = ro(t, r) gives the value of radial coordinate at time t = 0. Using 



dr vir) 
dr — vat = — — aro = — — - aro , 
oro v{ro) 

the line element in preferred (/* = 0) and comoving {g^i = 0) coordinates (t, vq) reads 



(54) 



drl + r^dn'^ 



(55) 



where r = r(t, tq). Therefore, the evolution of a comoving volume centered along the tq- 
trajectory is given by 



v{r) 
v{ro) 



X sin 6 = hl^2D ^ ^ • 



(56) 



As in cosmology, it possesses a well-defined time dependence. In the present case, it is governed 
by the (shift) function v{r), and the FF trajectories r = r(t,ro). 

Because of spherical symmetry, our Gaussian action separates into actions at fixed angular 
momentum : St = S/^^ 5't(/, m), see eq. (ITTl) . Each of them contains two-dimensional fields 
4>i,mit,r) and '^i^rn(t,r). We here consider only the s-wave sector, and drop the 0,0 index. We 
work with the rescaled field = r0, \|/ = r and drop the tilde. Taking into account the 4D 
character of the problem and eqs. (IHTf [3H1) . the action reads 



Si 



dtdrr 



Ihll 



2 

g^Lv 



r 



II 



dtdr 
dtdr 



(( 
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A 
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vdr + drV , 



dk^) 



(57) 



By varying this action, there is no difficulty to get the equations of motion. One verifies in 
particular that the kinetic action of the (rescaled) \E''~ = J\.^2^'^ field is as in 2D fiat space when 
expressed in preferred coordinates t, tq: 



S^ = ^jdkjj dtdro[{dt^',;}' - {(kALvTiWk)' 



(5J 



This means that the Chelva field \E'^(t,r) is a collection of independent oscillators labelled by 
tq. In other words, \E'^(t, r) depends on t and r as 



(59) 



One also verifies l'^dfj,'i/l{t,r) = 9t^^(t; ro)|ro , because l'^d^rQ{t^r) = 0, by definition. These 
facts guarantee that eq. flH^ applies. 
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By a direct analysis of the equation of motion, one obtains two important properties. First, 
one verifies tliat tliere is no " dilution" of the dissipation rate as time passes. This is non trivial 
(see [33j for an a priori similar system in which a dilution does occur). This steady behavior 
is due to the fact that the Chelva fields form a 3 dimensional dense set (as opposed to the 
discrete lattice of [33]) which allows the various to cancel each other in the determination of 
the dissipation rate. Moreover one also verifies that stationarity is preserved, not only in the 
dissipative aspects, but also in the driven properties govern by the kernel N, eq. ( 166|) . This is 
non trivial either since he depends on t and appears in the action S'^r^ when working with the 
rescaled fields . 

With these two results, one can study the impact of dissipation on Hawking radiation by 
applying the technics of pTO] which were further developed in [31][S2] to study the impact of 
dispersion. When the two scales are well separated, i.e. when k/A^v -C 1, where n = drV 
evaluated at the horizon v = —1, one can work in the near horizon approximation, wherein 
V = —1 + Kx, and in the p-representation. In this representation, the velocity profile becomes 
a non-trivial operator {) = — 1 — indp, in fact the only one. When working at fixed Killing 
frequency ujk, one obtains the radial propagation of the interacting degrees of freedom in terms 
of an adiabatic evolution of Green functions and kernels written in the ^-representation. These 
are given in Appendix B, with the frequency ut = idt of that Appendix replaced by the FF 
frequency Q = idt + ivdr = ujk — vp. We are planning to report on this soon. 
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4 Appendix A: Two-point correlation functions 

Since the models we consider are Gaussian, the complex function of eq. (fT5l) governs all 
observables built with the Heisenberg field operator 0. To analyse it, as mentioned in the text, 
it is appropriate to study separately the commutator and the anti-commutator. 
We start with the simple part, the commutator 

G,it,t';p)S%p-p')^TT[pT[Mi),'PU^')]-]- (60) 

It possesses several properties. First, using eq. ( 12^ one sees that it decomposes into two terms, 
one due to the non-commuting character of 0*^, the other due to that of Second, since both 
commutators are c-numbers, it is independent of pt, the state of the system. Hence, for all 
Gaussian models, one has 

G,{t,t';p) = [<l)\t),<p^{t')]_ + [[dhdt2Gr{t,h)Gr{t',t2)D{t^,t2), (61) 
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where the " dissipative" kernel D{ti,t2) is given by 

D{h,t2) = gM g.ih) = ^^9^{tl)Gl,{t,,t2)g,{t2). (62) 

Notice how this kernel combines the various couplings constants and the non-commuting prop- 
erties of the environment. 

The third property is the most relevant for us (and perhaps also less often mentioned). To 
all orders in gi and for all values of gi, Qi (even with arbitrary time dependence), the following 
identity 

idtGcit,t';p)\t=t' = 1, (63) 

holds because it corresponds to the ETC of eq. (fTUjl . The 1 on the rhs is guaranteed by the 
Hamiltonian character of the evolution of the entire system -|- It is therefore this equation 
which replaces the constancy of the Wronskian that was relevant in the case of free evolution. 
Eq. ( I63l) is crucial for us because the operator (f)'^{t) exponentially decays in t — tin where is 
the moment when the interactions are turn on, since it is an homogeneous solution of eq. fl25|) . 
Hence the first term in eq. flUTj) decays as exp — 7(t + 1' — 2tj„). Therefore at late times with 
respect to tin in the units of the inverse decay rate 7~^, the non-commuting properties of field 
operator are entirely due to those of the environment degrees of freedom, In addition, 
they sum up exactly to 1 as shown in the above equation. 
We now analyze the anti-commutator, 

Ga{t,t';p)5'{p-p')^TT[pT{Mt)AUt')U]- (64) 

When the density matrix factorizes as px = Pifi before the interactions are turned on (as it 
is the case in the "free" vacuum), Ga also splits into two terms. 

Gait, t'; p) = Ti[p^{ct>\t) , (t>\t')}+ ] + jj dt,dt2 Grit, t,)Grit', t2)iV(ti, ^2) • (65) 

The first term depends only on the initial state of 0. Similarly, the driven term depends on the 
state of the environment through the "noise" kernel 

Nitut2) = SiS,Tr[p^ {(^i(ti)^°(ti) ,^7,(t2)^°(t2)}+] . (66) 

As for the commutator, in the presence of dissipation, the first term exponentially decays, 
expressing the progressive erasing on the information contained in the initial state of 0. At late 
times therefore, as one might have expected, it is the state of \1/ which fixes the anti-commutator 
of 0. This allows to remove the restriction that initially the density matrices factorizes. If one 
is interesting by late time behaviour, only matters. 

In brief, two important results are been obtained. First, even when reduces to its driven 
term, the second term of eq. fl24l) . it still behaves as a canonical degree of freedom in that 
it exactly obeys eq. (j63ll . Second, only two (real) c- number quantities determined by the 
environment govern the two-point functions of 0, namely the kernels D and A^ of eqs. (|62l) and 
( l66l) . Therefore the set of environements (Gaussian or not) possessing the same kernels will 
give rise to the same 2pt functions for 0. Hence they should be viewed as forming an equivalent 
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class. The degeneracy can be lifted by considering correlations with observables containing the 
operators or higher order correlations functions of (j) (for non-Gaussian environements), two 
possibilities we shall not discuss in this paper. 

In order to be able to compute Gc and two different routes can be adopted. When gi{p) 
and VL{p) are constant, one should work in Fourier transform because the equations can be 
algebraically solved, in full generality. Instead when gi{p) and/or VL{p) are time-dependent, as 
it will be the case in expanding universes, it becomes imperative to choose the set of \l/i's and 
their frequency ^^{p) so as to simplify the time dependence of the equations. In the body of 
the article, we proceed with time dependent approach since it will allow us to combine general 
covariance and dissipation. We present below the Fourier analysis which is straightforward 
and well known [221 [2S] • We encourage the reader unfamiliar with the treatment of dissipative 
effects in Quantum Mechanics to read it. 

5 Appendix B : 

Stationary states, and vacuum 2-point functions 

In this Appendix, we present the relationships between Gc-, Ga which always hold in stationary 
states. In these cases, the Green functions G^ Ga and the kernels N are functions of 
t — t', and are related to each other in a fundamental way, generally refered as a Fluctuation- 
Dissipation relation. We briefly explain its origin and its physical implications in the present 
context. We start the analysis by with the most basic object: the retarded Green function Gr- 

5.1 The retarded Green function 

The Fourier transform of eq. ([21]) gives 

i-u' + ul)<i)p{u) = j^Mp^ip,^) + ^i9Kp)R°i^;p)M^) > (67) 

where 

R1{u;;p)={-{u; + ^er + n^{p)))-\ (68) 

is the Fourier transform (defined as in eq. f[TUl) ) of the retarded Green function of \I'i. As usual, 
its retarded character is enforced by the imaginary prescription of the two poles to lay in the 
lower half plane (e > 0). The solution of the above equation is 

0p(^) = 0j(^) + Gr{uj,p)J:,g,{p)^nP,^) , (69) 

where the Fourier transform of the retarded Green function of 0, the solution of eq. ( l25ll . always 
takes the form 

Griuj^p) = {-{u + lef + ujI + S,(cu,p))"' . (70) 

All effects of the coupling to the \l/j's are thus encoded in the (retarded) self-energy Y^riyj^p). 
For Gaussian theories, it is algebraically given by 

E,(^,p) = -S,^7.'(p)i?°(^,p). (71) 
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The dissipative effects are governed by the imaginary part of Ilr{uj,p)- In the present case, one 
has 

2ImS,(^,p) = -J:,g^ip)Gl,iuj,p) = -D{uj,p) . (72) 

To get the first equahty we have used the fact that in stationary states the retarded Green 
function and the commutator are related by 2\m.Gr{uj) = Gc{uj) for all degrees of freedom, free 
or interacting. In the second equality, we have introduced D{uj), the Fourier transform of the 
kernel of eq. (l62ll . 

Several observations should be made here. First, from eq. (l68ll . we obtain that D{uj) is 
proportional to T,igf6{uj — Qi). Therefore there is no dissipation for lower frequencies than the 
lowest value of fli. This simply follows from energy conservation. Second, to obtain "true" 
dissipation, D{u,p) should be a continuous function and not a sum of delta. This can only 
happen when the form a dense ensemble. In the body of the paper, we shall thus replace the 
discrete sum on i by an integral on a continuous variable, k. We shall not consider the discrete 
cases even though these could display interesting properties. Third, from a phenomenological 
point of view, only D{u!,p) matters. Hence we cannot separately know what is the spectrum of 
the environment, which is given by R°{uj,p), and what is the coupling strenght gf{p)- This is a 
good thing, because when working in time-dependent settings, we shall exploit this equivalence 
by choosing the simplest model of \E'i's which delivers the kernel D{uj,p) we want. 

It is also worth noticing that the dispersive (real) effects are not directly related to D (or 
A^). These are governed by the even part of Sj,(co',p) which is given by 

ReS, = /- 73 

J Ztt uj — uj' 

where the integral should be understood as a principal value. This integral relationship explains 
why one often founds that dispersive effects appear before dissipative effects (for increasing ti;). 
We also learn that the dispersive models studied in the litterature are incompatible with the 
above relations since they assume both ReS^ 7^ and ImS^ = 0. Therefore these models 
cannot be viewed as resulting from dynamical processes. 



5.2 Fluctuation-Dissipation relations and vacuum self-energy 

In this subsection, we derive the relationships between Gc, Ga and J^p which exist in the true 
(interacting) ground state. 

In interacting theories, the only stationary states are thermal states, see e.g. [32]. In these 
states, the Fourier transform of D and are related by 

N{iu) = D{uj) coth(/?cu/2) , 

= D{u)sign{uj)[2n{\u\) + l]. (74) 

In the second line, n{uj) is the Planck distribution. It gives the mean occupation number of 
quanta as a function of the frequency (measured in the rest frame of the bath). The above 
relation directly follows from the fact that the individual commutators and anti-commutators 
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of the free fields obey this relation, as any free oscillator would do. It implies that the 
Fourier transform of Gc and Ga are also related by 

Ga{uj) = G,{u) sign(cu) [2n{u) + 1] . (75) 

It should be stressed that this equation is exact, i.e. non-perturbative, and valid for all theories, 
Gaussian or not. (It indeed directly follows from the cyclic properties of the trace defining 
Gpit,t') =Tr[e-^^-<^(t)</.(t')]). 

For Gaussian models, there exists an alternative direct verification of eq. fl75|l . In fact, it 
suffices to note that in steady states the decaying terms of eqs. flUTl) and flU3|) play no role, and 
that the Fourier transform of the driven terms are respectively given by 

G,{u) = \Gr{u;)\^ Diu) , (76) 
Gaioo) = \Gr{uj)\^N{uj), (77) 

since the Fourier transform of the retarded Green function obeys Gr{uj) = G*{—uj), see eq. 
(169|) . Irrespectively of the complexity of Gr, i.e. irrespectively of the functions gi{p), ^i{p) and 
the set of the \Efj fields, Gc and Ga are thus related to each other by the FD relation (!73|) . 

These universal relations will be relevant for to infiationary models wherein only the ground 
state contributes. In particular, they imply that in the true vacuum, i.e. when n{u) = 0, Gc 
and Ga are exactly related by Ga{uj) = Gc{uj) sign(co'). Hence the Wightman function 

Gw = ^{Gc + Ga)=9{u;)Gc, (78) 

is determined by the commutator and contains only positive frequency, as in the free vacuum. 
Equations ( 1761 1771) also allow to compute the vacuum self-energy of the Feynman Green function 
in eq. ([2]). For Gaussian models it is given by 

2ImEir(tu) = -D{u) sign(cj) = 2Re E^^f G'^^^(cj) , (79) 

where the p dependence has not been not explicitized, and where we have taken the real part 
because of the i in the numerator of eq. ([2]). With the last equality we recover the fact that in 
the vacuum, it is sufficient to consider Feynman Green functions. In non-vacuum states, and 
in non-stationary situations, this is no longer true, thereby justifying the use of the in — in 
machinery (also called Schwinger-Keldish formalism), see e.g. [27] . 

Before specializing to a specific class of models giving rise to dissipation at high frequency, 
we make a pause by asking the following important question: What should be known about 
the \E'j fields to get eqs. fl76| [771 178|) ? We have proven that it is sufficient for the \E'j's to be 
canonical fields, but is it necessary ? 

The answer is two fold. On one hand, the \t'j cannot be stochastically fiuctuating (i.e. 
commuting) quantities because this would lead to a violation of eq. ( 17il) that would imply the 
violation of eq. (JTSl) and the ETC eq. They cannot be either a combination of quantum 

^ This constitutes the simplest proof that it is inconsistent to couple quantum variables to stochastic (or 
classical) ones. If one does so, the ETC of the dressed quantum variables will always be dissipated after a time 
of the order of 7"^. One can therefore view the experimental evidences for the ETC of some degrees of freedom 
as a very strong indication that all dynamical variables in our world are quantum mechanical in nature. This 
line of thought has been used by W. Unruh to prove that gravitational waves must be quantized. 
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and stochastic quantities because this would still lead to a violation of the ETC. Hence they 
must be built only from quantum (canonical) degrees of freedom. 

On the other hand, the \E'j's can be composite operators]^ i.e. polynomials of some (un- 
known) canonical fields. Indeed, their commutators would still be all related to their anti- 
commutators by the FD relation eq. fl7^ . and this even though both depend non-linearly 
on n{uj) in non-vacuum states. The difference with Gaussian models is that these non-linear 
operators possess non-vanishing connected higher order correlation functions. Hence, the self- 
energies Sr, Si? will contain higher a series in powers of gf, and not just a single quadratic 
term as in eq. (ITTl [79]) . Nevertheless these higher loops corrections preserve the validity of eq. 
(!75|) in the ground state, as well as that of eqs. (17^ 1771) in any thermal state, when properly 
undertood, i.e. with D now defined by -2ImSr (as the effective dissipation kernel), and N 
related to it by the FD relation. 

In brief, we have reached/recalled the following results. Firstly, the q-number combination 
S = T.iQi'^i, the fluctuating source term of 0, must obey the FD relation fTMl) . This can either 
be postulated, or better, be viewed as resulting from the fact that S is entirely made out of 
quantum degrees of freedom. Secondly, to lowest order term in g, the self-energy can be obtained 
by treating S as a quantum Gaussian variable, whatever its composition may be. Thirdly, when 
dealing with non-Gaussian theories, once having computed T,r{u!), the resulting equations for 
the 2-point functions have the same structure and the same meaning as in Gaussian theories, 
with D replaced by -2ImSr. Therefore, the entire phenomenology of two-point functions can 
be described with Gaussian settings. 

5.3 The double limit: g'^T oo followed by g'^ 0. 

To perform a phenomenological analysis, we need to understand how the theory behaves in 
transitory regime from dissipative to free propagation. Similarly, to study primordial spectra 
in inflation or Hawking radiation, we also need to understand how free motion emerges as the 
proper frequency get red-shifted. It is therefore relevant to study the behavior of the two-point 
function in the following double limit. 

One first takes g^T oo, where T = t — tin, tin being the moment when the interactions 
are turn on, and t the moment when one studies the field properties. This limit implies that 
the decaying term in eq. 0241) plays no role. Therefore, near time t, the Heisenberg field 0(t) is 
a composite operator which only acts in the Hilbert space of S. 

Secondly, one considers the "free" limit ^ of that composite operator. One could 
naively conclude that Gc and Ga of eqs. (1761 1771) would vanish since both D and are 
proportional to g^. However, this is not the case, because the prefactor in these equations, 
IGrp, is singular in this second limit. In fact, one verifies that it scales in l/g"^ in such a way 
that, in the (interacting) vacuum, one always recovers 

Gw{uj)g2^o = TT—'^'^^i^ - ^p) ■ (80) 
Zujp 

This is the standard vacuum fluctuations of a free oscillator of frequency ujp. 
^'^l am grateful to Albert Roura for bringing my attention to this question. 
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Two important lessons have been reached. First we learned is that even though (j) acts only 
on the S-Hilbert space, when g"^ —>■ 0, it behaves as if it were a free mode possessing its own 
Hilbert space, with no reference to S-dynamics. Secondly, in spite of this, the quantum state 
in the would be Hilbert space is still exactly that of S. Therefore, in stationary situations, 
the only "souvenir" kept by the composite operator is the equilibrium distribution n{uj) 

inherited from its parents. 

Let us now emphasize that the above limit is relevant for non-Gaussian models as well. 
Indeed, in the limit ^ 0, there will always be a value of g"^ sufficiently small that the model 
can be well approximated by a Gaussian model. Therefore the behavior of the 2-point functions 
in the transitory regime from dissipation to free propagation can be analyzed without restriction 
by studying Gaussian models (at least in the quasi-static limit). 



6 Appendix C: 

Dissipative effects above Aiv- The Phenomenology 

We now have all the tools to understand models giving rise to dissipation in the vacuum above 
a critical energy scale Alv- 

If one considers only stationary situations (i.e. static metrics and stationary states), and 
if one adopts a phenomenological point of view, one can simply choose the function D{uj,p) 
entering eq. 0721) and eq. (!7H|) as one wishes. There is indeed no restriction on D{uj,p) besides 
its constitutive properties, namely being odd in u and giving rise to poles in Gr all localized in 
the lower half u plane. In this we have reached our first aim, namely identify how to generalyze 
the free settings so as to incorporate some arbitrary dissipative effects. 

We can thus consider the dispersive models which correspond to those defined by eq. 
They are charaterized by a single term giving rise to dissipation above Alv Iii the vacuum, 
they are fully specified by the imaginary part of the (retarded) self-energy 

-lmJ:^-)(u;,p) = -^p'(-^) =2^7n. (81) 
In these models, the decay rate (inverse life time) on the mass shell is 

/ \ 271+1 

To verify it, assuming that ReS^ = 0, the two poles of Gr{uj) in eq. (!70l) are located in 

u±{p) = ±^ul-^^-i^. (83) 

From this, by inverse Fourier transform G^{uj), one obtains that the decay rate is indeed 7 in 
the underdamped regime, for uj"^ > 7^- In the overdamped regime, for 7^ > i^p, the decay rates 
of the two independent solutions of G^Vd = are F-t = 7 ± 1^7^ — oj^. 
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One thus obtains the following behavior as p grows. For p <C A^y, uj± ~ p, and one has a 
free propagation which is slightly damped with a life time in the units of the frequency given 
by (Aiy/u;)"'+^ ^ 1. Instead the opposite regime of high momenta p ^ A^y, deep in the 
overdamped regime, the two roots lj± are real and the notion of propagation (in space-time) is 
completely absent. In anticipation to what will occur in inflation or in black hole physics, we 
invite the reader to study the migration of the poles of Gr when extrapoling backwards in time 
a mode, i.e. as p increases. (Remember that the physical momentum of a mode in cosmology 
is p{t) = po/a{t) where po is the norm of the conserved comoving wave vector, whereas near 
a black hole horizon one finds p{r) = u/x where x = r — is the proper distance from the 
horizon measured in a freely falling frame, and u the conserved Killing frequency measured 
asymptotically.) 

One could of course generalize the above class by considering in eq. (IHTl) polynomials in 
p dimensionalized by different UV scales. However, unless fine tuning, the phenomenology of 
the transition from the IR dissipation-free sector to the dissipative sector will be dominated a 
single term. One should also consider the possibility that ImS strictly vanishes below a certain 
frequency ^i. This would be the case when the spectrum of the fields possesses such a gap, 
see the remarks after eq. ( 1721) . 

Having the phenomenology of dissipative unitary models under control (with dispersive and 
dissipative related by Kramers relations, see eq. ( 173|) ) one could confront particle and astro- 
particle physics data and put lower bounds on A^y for each n, in analogy with what was done 
for (pure) dispersion in 0. 
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